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A.  Background 

Under  this  Phase  I  SBIR  effort,  Mission  Research  Corporation  (MRC)  and  the  University  of 
Michigan  (UM)  undertook  a  feasibility  study  involving  the  application  of  recent  innovations  in  the 
finite  element  (FE)  method  to  the  problem  of  modeling  conformal  antennas  that  utilize  ferrite 
materials.  This  final  technical  report  presents  a  comprehensive  technical  summary  of  the  work 
undertaken  during  Phase  I. 

Ferrite  materials  have  been  suggested  by  members  of  the  antenna  community  as  a  unique  antenna 
loading  material  primarily  due  to  its  tunable  electrical  properties.  In  particular,  by  appropriate 
application  of  a  static  magnetic  bias  field,  a  ferrite  sample  can  exhibit  a  wide  variety  of  electrical 
properties  from  nearly  transparent  to  nearly  opaque  without  any  mechanical  intervention.  Ferrites  can 
also  be  magnetically  tuned  to  preferentially  interact  with  a  given  electromagnetic  field  polarization. 
For  example,  the  sample  can  appear  to  be  opaque  to  one  polarization  while  transparent  to  the 
orthogonal  polarization.  Such  properties  can  be  exploited  by  antenna  and  low  observable  (LO) 
designers.  In  this  report,  preliminary  research  concerning  such  properties  will  be  presented. 

One  difficulty  in  designing  such  antennas  involves  the  dearth  of  appropriate  analysis  and  design  tools. 
Traditional  analysis  methods,  utilizing  the  popular  method  of  moments  (MoM)  procedure  (i.e. 
integral  equations),  have  proven  unsatisfactory,  principally  due  to  the  fact  that  the  static  magnetic 
field  cannot  be  properly  represented  using  typical  MoM  techniques.  Practical  biasing  magnets 
produce  a  non-uniform  static  magnetic  field  and  hence  the  ferrite  sample,  though  physically 
homogeneous,  appears  to  be  electrically  inhomogeneous.  Volume  MoM  analysis  programs,  which 
can  indeed  model  such  electrical  inhomogeneities,  are  inefficient  and  difficult  to  formulate  and 
program  for  anisotropic  materials  such  as  ferrites. 

The  finite  element  method  is  uniquely  suitable  for  analyzing  the  properties  of  ferrite  materials.  This 
is  due  to  the  fact  that  the  finite  element  method  is  typically  formulated  in  the  frequency  domain,  and 
hence  can  readily  accommodate  the  frequency-dependent  properties  of  ferrites,  hi  addition,  since  the 
finite  element  method  requires  no  Green’s  function  for  the  ferrite  volume,  its  formulation  and 
implementation  is  particularly  easy.  Also,  since  the  finite  element  method  is  based  on  partial 
differential  equations  (PDE),  it  possesses  a  locality  property  in  that  only  field  expansions  in  the  near 
vicinity  of  one  another  couple  directly.  All  of  these  features  lead  to  the  conclusion  that  the  finite 
element  method  is  well-suited  for  ferrite  material  modeling. 

The  focus  of  this  Phase  I  SBIR  topic  involves  the  investigation  as  to  whether  the  FE  method  can  be 
used  to  model  ferrite-loaded  antennas.  In  addition,  it  seeks  to  begin  to  answer  the  question,  “What  are 
the  advantages  and  disadvantages  of  ferrite-loaded  antenna?”  Note  that  in  this  report,  a  ferrite-loaded 
antenna  does  not  include  the  ferrite  bead  loads  commonly  used  to  alter  the  performance  of  certain 
antennas.  Rather,  the  ferrite  load  is  either  a  substrate  or  superstrate  material. 

In  this  report,  the  finite  element  method  will  be  shown  to  be  well-suited  for  ferrite  antenna  analysis. 
Also,  novel  uses  of  such  antennas  will  be  explored  including: 
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1 .  Frequency  agile  antennas. 

2.  Polarization  agile  antennas. 

3 .  Radar  cross  section  (RCS)  reduction  techniques. 

We  begin  this  report  with  the  finite  element  formulation  for  ferrite  materials  followed  by  a  discussion 
of  MRC’s  and  UM’s  efforts  during  Phase  I. 

B.  Finite  Element  Formulation  for  Anisotropic  Materials 

The  finite  element  formulation  presented  herein  is  very  general.  It  does  not  assume  any  particular 
element  shape,  expansion  (or  shape)  function,  or  mesh  termination  condition.  For  practical  analysis, 
the  formulation  will  need  to  be  specialized.  Since  during  Phase  I,  both  MRC  and  DM  utilized  brick 
elements  and  a  boundary  integral  (BI)  mesh  closure,  these  special  cases  will  be  presented  following 
the  general  discussion.  Note  that  both  MRC  and  UM  chose  to  use  bricks  and  a  BI  termination  to 
leverage  previous  work  and  thus  be  able  to  quickly  investigate  the  special  properties,  features,  and 
challenges  associated  with  ferrite  materials.  A  more  general  treatment  is  planned  for  Phase  II. 

B.1  General  Finite  Element  Formulation  for  Anisotropic  Materials 

The  finite  element  formulation  begins  with  Maxwell’s  equations 

Vxe*‘=— jk^ir*  (i) 

V  x  H”*  =  +jk0Z0erE“t  (2) 

where  the  interior  (e.g.  within  the  computational  volume)  electric  and  magnetic  fields  are  denoted  by 

Emt  and  Hmt,  respectively.  The  relative  permittivity  (er)and  permeability  (pr)are  tensors  and  can 
vary  on  an  element-by-element  basis  within  the  computational  volume.  The  free-space  wavenumber 
and  impedance  are  denoted  by  k0  =  co^/p  e0  and  Z0^p0/e0 ,  respectively,  where  the  free-space 
permittivity  and  permeability  are  represented  by  e  0  and  p0 ,  respectively.  Pre-multiplying  (1)  by  the 

inverse  permeability  tensor  (pr_1),  taking  the  curl  of  (1)  and  utilizing  (2),  the  second-order  PDE  is 
obtained 


V  x  [p;1  •  V  x  E”‘ ]-  k*er  -Emt  =  -[jk0Z0 J“p  +  V  x (p;1  •Mimp)]  (3) 

This  is  the  vector  wave  equation  in  toms  of  the  total  electric  field  within  the  computational  volume, 

V  which  is  illustrated  in  Figure  B-l .  In  (3),  J imp  and  M“np  are  the  impressed  electric  and  magnetic 
current  within  the  computational  volume,  respectively. 
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Figure  B-l .  Illustration  of  the  computational  volume  (V),  bounding  surface  (S),  and  surface  at  which 
the  radiation  condition  is  applied  (so) . 

Throughout  this  report,  total  electric  fields  will  be  the  unknowns  since  the  BI  mesh  closure  used 
during  Phase  I  involved  a  half  space  Green’s  function  (e.g.  cavities  recessed  in  an  infinite  metallic 
plane)  and  hence  a  total  magnetic  field  formulation  is  impractical.  This  is  not  to  say  that  a  magnetic 
field  formulation  will  not  be  preferable  for  ferrite  materials.  Indeed,  the  fact  that  a  magnetic  field 

formulation  precludes  the  need  for  an  inverse  permeability  tensor  (p  )  is  desirable  due  to 
singularities  in  that  tensor  near  the  forced  precession  frequency.  Hence,  during  Phase  II  it  will  be 
important  to  examine  the  benefits  of  a  magnetic  field  formulation  utilizing  an  artificial  mesh 
termination  condition  such  as  perfectly  matched  layers  (PML) 

Two  methods  are  typically  used  to  convert  the  continuous  domain  wave  equation  represented  by  (3) 
and  a  discrete  approximation  that  can  be  solved  on  a  digital  computer,  namely:  the  method  of 
functionals  and  method  of  weighted  residuals.  Throughout  this  report,  the  method  of  weighted 
residuals  will  be  used  principally  due  to  the  fact  that  it  is  familiar  to  electrical  engineers  who  have 
used  the  method  of  moments  before.  The  method  of  weighted  residuals  begins  by  multiplying  (3) 

with  an,  as  of  yet  arbitrary,  weight  function  (w*)  that  has  support  only  over  one  finite  element  (e). 
There  are  N  such  unique  weights  corresponding  to  the  number  of  unknowns  in  the  final  system. 
After  multiplication,  (3)  is  integrated  over  the  entire  computational  volume  (though  of  course  the 
finite  support  of  the  weight  functions  will  limit  the  integration  to  the  element  associated  with  the 
weight  function)  to  get 

|v  {wie  • V  x [pr1  •  V  x  Emt ] -  k'W'  •  Ir  •  Emt }  dV  = 

-fvWT  •[jkoZoJ”P  +Vx(pr-'  •  Mimp )] dV  =  f®’”*  (4> 

where  the  right-hand  side  (interior  excitation)  is  represented  by  f',mt  for  the  remainder  of  this  report. 
This  equation  is  not  convenient  for  solution  since  the  electric  field  receives  two  derivatives  and  the 
resulting  formulation  is  unsymmetric.  Therefore,  it  is  advantageous  to  transfer  one  derivative  from 
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the  unknown  electric  field  to  the  weight  function  in  the  first  integrand  using  the  first  vector  Green’s 
theorem 

£w;  •  V  x  •  V  x  E"  ]  dV = Jv  ( V  x  W,' )  ■  (5;'  •  V  x  E“  )  dV  -  £n  ■  (w,‘  x  p;'  ■  V  x  E“  ])  dS 

(5) 

Using  (1)  and  a  vector  triple  product,  (5)  becomes 

jw”  •  V  X  p;'  •  V  x  Em]  dV  =  fv(v  x  W*)  •  p;1  •  V  X  Eta)  dV  -  jkoZo^W,8  •  (n  X  Hto)  dS  (6) 

In  this,  the  surface  term  relates  die  tangential  electric  field  (\V*)with  the  tangential  magnetic  field 
(n  x  Hm)on  the  surface  enclosing  the  volume  (V)  where  n  is  the  outward  directed  normal  to  the 
surface  S  (dV) .  Using  (6)  in  (4),  the  weak-form  of  the  vector  wave  equation  is  obtained 

Jv[(v  X  w ,•).  jT;1  -(v  X E"‘ ) - k*W*  •  8r  -Emt ]  dV- jk0Z0  j>W,e  -(n  x  Hmt )  dS  =  f*’“‘  (7) 

This  equation  satisfies  (3)  in  the  weak  or  integrated  sense  rather  than  at  every  point  within  the 
computational  volume. 

Finite  element  programs  are  typically  divided  into  two  groups  based  upon  how  the  surface  magnetic 
field  is  related  to  the  surface  electric  field: 


1 .  locally  through  a  approximate  condition  such  as  absorbing  boundary  conditions  (ABC)  or 
perfectly  matched  layers  (PML),  or 

2.  globally  via  a  boundary  integral  involving  a  dyadic  Green’s  function  for  the  external 
region. 


During  Phase  I,  both  MRC  and  UM  used  the  later  approach  and  its  formulation  will  be  presented 
next. 


B.2  Finite  Element-Boundary  Integral  Formulation  for  Anisotropic  Materials 

The  finite  element-boundary  integral  formulation  begins  by  postulating  a  magnetic  field  integral 
equation  (MFIE)  for  the  exterior  region  such  as 


x^V  x  G-n' x  H^fr')! 
nxjH^rJ+H^r)] 


dS,-jk0Y0|snxGn'xEex‘(r')]dS  = 


(8) 


where  the  field  quantities  are  now  in  the  exterior  region  (e.g.  between  S  and  S0in  Figure  B-l),  the 
prime  variable  (r')  refer  to  the  source  point  while  the  unprimed  variable  (r)  refer  to  the  test  point, 
and  Y  =  1/Z0 .  The  excitation  (the  right-hand  side)  is  separated  in  term  of  an  incident  magnetic 
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field  (Hmc)and  a  reflected  magnetic  field  (h^ ) .  For  geometries  involving  an  infinite  metallic 
structure,  this  field  representation  is  used.  For  geometries  involving  finite  structures,  the  reflected 
field  is  omitted  and  for  radiation  analysis,  both  the  incident  and  the  reflected  fields  are  omitted. 

This  MFIE  is  then  tested  with  a  separate  vector  test  function  (q*  ) 

^Q'  •nxH“t  dS-££ Q'  •  n x [v x S- n' x H“‘ ] dS'dS- 

jk0Y0 ££,  Q*  •  n  x  p-  n'  x  E”1  ]  dS 'dS  =  £Q'  •  n  x  [H"0  +  ]  dS  =  f*'m 

and  the  interior  and  exterior  fields  are  matched  across  the  interface  using  either  point-by-point 
coupling 


nx[Eint-Eext]  =  0  and  n xIh^ -H“]  =  0  (10) 


or  in  an  average  sense 

j>gQ*  •  n x [e”1*  - E°“]  dS  =  0  and  j>Q‘ -hx^ -H“]  dS  =  0  (11) 


Although  the  formulation  is  appropriate  for  general  structures,  such  as  the  one  illustrated  in  Figure  B- 
1,  for  Phase  I  both  MRC  and  UM  used  a  specialized  FE-BI  formulation  for  cavities  recessed  in  an 
infinite  ground  plane,  such  as  the  one  shown  in  Figure  B-2,  to  save  on  computer  resources. 


Figure  B-2.  Illustration  of  a  cavity-backed  aperture  located  in  an  infinite  ground  plane. 

This  is  due  to  the  fact  that  a  particularly  efficient  FE-BI  computer  program  was  developed  by  UM  to 
analyze  such  structures  and  using  this  program,  both  MRC  and  UM  were  able  to  quickly  investigate  a 
large  number  of  ferrite  antenna  configurations.  Such  efficiency  is  important  for  performing  an 
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exhaustive  investigation  of  ferrite  antenna  capabilities,  idiosyncrasies,  and  potential.  This  special 
formulation  allowed  the  use  of 

zxH6*  =zx[Hinc+Hre]  +  jk0Y0£zxSxz-E“t  dS'  (12) 

where  we  have  assumed  that  the  metallic  plane  is  at  z=0  and  hence  the  normal  unit  vectors  become 
n  =  h'  =  z.  If  identical  expansion  functions  are  used  for  both  the  exterior  surface  mesh  and  the 
interior  volume  mesh  (that  is  the  volume  shape  function  reduces  to  the  surface  shape  function  as  the 
evaluation  point  approaches  the  surface)  which  results  in  an  implicit  enforcement  of  (10),  then  (12) 
can  be  substituted  into  (7).  This  results  in  the  FE-BI  equation  used  in  Phase  I 

fv [(V  x  W* ) •  f*r“‘  -(v  xEmt)-koW'  er  •Emt]dV+  k2  j>  £  W*  •  [z x  G  x  z] •  Emt  =  f'-”*  +  f *■"“ (13 

) 


For  cavities  recessed  in  a  metallic  plane,  the  half-space  dyadic  Green’s  function  should  be  used  in 
(13) 


G  =  2G0 


=  2 


’=  vv’ 

’=  vv’ 

I+ p" 

L  J 

II 

o 

O 

I+k2 

L  Ko  J 

ejk0|r-r'| 

27t|r  —  r'| 


(14) 


The  FE-BI  equation  (13)  is  complete  and  it  only  remains  to  specify  the  finite  elements  and  then- 
associated  shape  functions. 


B.3  Finite  Element  Method  using  Brick  Elements  for  Anisotropic  Materials 

The  finite  element  method  can  utilize  many  different  finite  element  shapes,  several  of  which  are 
shown  in  Figure  B-3 . 


Right  Angled 
Brick 


Tetrahedron 


Skewed  Brick 


Curvilinear  Brick 


Cylindrical  Shell 


Figure  B-3.  Various  finite  element  shapes. 


Mission  Research  Proprietary  Information 


Mission  Research  Proprietary  Information 


10 


For  Phase  I,  both  MRC  and  UM  used  brick  elements  such  as  the  one  shown  in  Figure  B-4. 


8  7 


Figure  B-4.  Brick  finite  element  illustrating  the  six  nodes  and  twelve  edges  that  form  the  brick. 

These  elements  have  two  advantages : 

1.  the  finite  element  volume  element  matrices  can  be  solved  in  closed  form  and  in 
fact  only  a  small  number  of  these  interactions  need  be  computed  and  stored  for 
each  layer  regardless  of  the  material  variations. 

2.  they  result  in  a  uniform  discretization  of  the  surface  aperture  and  this  results  in 
boundary  integral  matrices  that  are  Toeplitz. 

Each  of  these  items  are  discussed  separately. 

B.3.1  Anisotropic  Brick  Element  Matrices 

The  element  matrices  for  bricks  and  anisotropic  media  can  be  written  as 

A">  =  Jv (v x w, •)•£■'  (v x w;) dV,  (15) 

A'2,=Jv.W'-?r-W*dVe  (16) 

where  the  inverse  permeability  and  permittivity  tensors  are  given  by 
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The  interactions  associated  with  (15)  are  given  by 


KK~  KK  -  „  K-  „  1-  „T 

6h*  ^'zz^1  +  av.®  ^yy^2  +  a  l*zy^3  +  A  ^yz®^ 


6h! 


A"1  =  +^-(I„k4  +^P„k6 

a™  =  -  yM„Kj  +^|rP„Kj  +^-ff„K6+^-it„Ks 


h*  h;  h'  h:h; 

A<”  =  -jP^K,  +-f  ii^K.  +-j-S„K6  +^fP„K 


fcyx 


T 

xyA^3 


KK~ 


KK 


K 


*2  =  6h^Sl“Kl  +  tf^*2  +  4  f +  4 


...  h'  _  hjh'  _  h*  h’ 

A">  =  — |-n„Ks +-^ri*,Ks  +-^n„K4 +^-h„k, 


...  h!  _  h!h:  _  h*  h* 

A">  =  -f  |l„Kj  +-^A|l„K3  +^-(*wK4  +~HyIK6 


h'  ^  h'h*  _  h'  „  h;  _ 
Aj>=— ^M„KJ+-^r|lwK,+-^n„K.+-^|l,K1 


X 

.  Ct  e 


h;h;_  h'h'  h-  h*  T 

A“  +"sh^,1“Kj+^>i>"K3+^^K’ 

x  y 


(19) 


while  those  associated  with  (18)  are  given  by 


a(2)  =  h:h;hz 


8 

— —  Kn 

!s.L 

24  Ll 

gxz  j  T 

24  Ll 

36  0 

8 

8 

e 

—  L* 

—  K0 

— L, 

24  1 

36  0 

24  1 

—  L 

24  Ll 

f-L 

24  L' 

—  K0 

36  0 

(20) 


where  the  superscript  “T”  denotes  a  matrix  transposition.  The  brick  edge  lengths,  |h®  ,h* ,  h*  J  are 
illustrated  in  Figure  B-4.  The  component  matrices  used  in  (19)  and  (20)  are  given  by 


Mission  Research  Proprietary  Information 


Mission  Research  Proprietary  Information 


"4 

2 
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1] 

2 

4 

1 

2 

K,  = 

2 

1 

4 

2  1 

1 

2 

2 

4J 

2  1-2-1 
1  2-1-2 
K2=  _2  -1  2  1 
-1  -2  1  2 
"1-11  -1" 
-11-11 
K4=  1  -1  1  -1 
-11-11 
= 1  1-1  -f 

1  1-1-1 
K--l  -1  1  1 

-1-11  1 


K3  = 


2 

-2 

1 

-1 

-2 

2 

-1 

1 

1 

-1 

2 

-2 

-1 

1 

-2 

2 

-1 

-1 

1 

1  " 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

1 

1 

-1 

-1 

Ks  = 


2 

1 

-2 

-2 

-1 

2 

1 

2 

-1 

-1 

-2 

1 

"2 

1 

2 

2 

1 

2 

1 

2 

1 

1 

2 

1 

B.3.2  Boundary  Integral  Matrix  Entries 

The  boundary  integral  in  (13)  can  be  represented  as 

IBI  =  2|Jg;  [wte  •  (z  x  cf0  x  zj  Wse ]  dS'dS  =  IBI(,)  +  IBI(2) 
where  G0  is  the  free-space  dyadic  Green’s  function.  From  this,  we  recognize  that 
IBI(,)e  =  -2jjs,  [wte  (z  X  X  z).  w;  ]  dS'dS  (23) 


lBI(2) ="^U{w‘{£xVV-^xiJ-w'Jds'ds  <24> 

where  R  =  |r-r'|  =  -J(x-x')2  +(y-y')2  +(z-z')2  .  In  Phase  I,  edge-based  brick  elements  were 
used  for  the  expansion  and  so  these  basis  functions  can  be  represented  by 
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w*  =T ^(y~y) 

(25) 

w;=^(x-,) 

as  the  evaluation  point  approaches  the  surface  (x, y,  z  h-»  o) .  In  (25),  the  parameters  x,  y , and  s 
are  chosen  based  upon  the  particular  local  edge  involved  in  the  calculation.  Using  (25)  in  (23),  the 
first  class  of  BI  interactions  become 


T  BI(1) 

Axx 


T  BI(1) 

yy 


yUiyU; 


vt  „t  „* 

yi  xi  yi  xi 


/ill  (y-y,)(y'-yj- 


-jkoR 


R 


y»x»yixi 


y|  xi  yf  xf 


JJJJ  (x-X,)(x'-x.): 


-jkoR 


dx'dy'dxdy 


dx'dy'dxdy 


(26) 

(27) 


where  the  limit  subscripts  “1”  and  “u”  refer  to  the  lower  and  upper  limits  of  the  integration  (e.g.  the 
boundaries  of  the  element). 


The  boundary  integrals  in  (26)  and  (27)  cannot  be  evaluated  in  closed  form  (as  was  the  case  with  the 
finite  element  contributions  given  above).  Rather,  these  integrals  must  be  evaluated  numerically  and 
special  care  must  be  taken  for  the  self  cell  (e.g.  when  the  source  and  test  point  can  coalesce).  In  these 
circumstances,  the  self-cell  can  be  subdivided  into  sub-cells  over  which  die  basis  function  is  assumed 
to  be  constant  and  closed  form  evaluations  for  the  resulting  sub-cell  integrals 


y.x,y;x‘. 


y,  x,  y,  x,  L 


R 


dx'dy'dxdy 


(28) 


are  available.  For  the  case  of  interactions  between  different  cells,  simple  numerical  integration  is 
sufficient. 


The  second  boundary  term  requires  further  manipulation.  Specifically,  applying  the  divergence 
theorem  twice  and  converting  one  of  the  gradient  operators  from  V  to  V'  (e.g.  V'G0  =  -VG0),  (24) 
becomes 


I“">  =‘iIrlX[v(J><w')v'(S!<w.) 


,-jkoR 


R 


dS'dS 


(29) 


Evaluating  (29),  using  the  basis  function  given  in  (25),  four  terms  are  obtained  based  upon  the  mix  of 
test  and  source  edge  orientation 
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jBI(2)  _ 


*t*k 


i“«>  =  + — iL 

2*j(h;h;) 

C’  =+ — t* — \ 

2ukS(h;hj) 


.jBI(c) 


jBI(c) 

jBI(c) 


tBI(2)  _  __  jjt  Ss 

27±2o(kY 


jBI(c) 


(30) 


Once  again,  the  IBI(C)  terms  in  (30)  must  be  evaluated  by  subdivision  and  closed  form  expressions  for 
the  self-cell,  and  via  numerical  integration  for  the  remainder  of  the  interactions. 

C.  Ferrite  Material  Model 

The  FE-BI  formulation  presented  in  Section  B  is  valid  for  arbitrary  anisotropic  materials.  The  focus 
of  this  SBIR  project  on  the  other  hand  is  a  anisotropic  material  with  special  properties,  ferrites.  These 
material  exhibit  dispersive  properties  that  can  be  manipulated  by  adjusting  the  strength  and 
orientation  of  a  static  magnetic  field.  That  is  the  material  “looks”  different  to  a  dynamic  field  as  the 
static  “control”  field  is  altered.  In  this  section,  the  ferrite  material  model  will  be  presented  along  with 
some  important  parameters  and  terms  that  will  be  used  throughout  this  report. 

The  most  common  representation  of  the  material  properties  of  a  ferrite  involve  a  scalar  permittivity 
(typically  fairly  large  ~9  to  16)  and  a  tensor  permeability  with  different  forms  depending  on  the 
orientation  of  the  applied  static  bias  field  (e.g.  this  field  is  polarized  along  the  direction  of  the  length 
of  the  magnet).  The  Polder  tensor  for  the  three  Cartesian  directions  are  given  by 


1 

0 

0 

Hr  = 

0 

H 

jK 

0  - 

jK 

H_ 

H 

0 

jK 

II 

ll=L 

0 

1 

0 

_-JK 
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o' 
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1 

x-bias 


y-bias  (31) 


z-bias 


where  the  p  and  k  terms  are  given  by 
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p  = 


1+2  2 
V  <»0-©V 


1+- 


fgfm 

fo-f2 


(0  to. 


K=l 


.©o  “©V 


ffm 

fn  -f2. 


Vi0 


(32) 

(33) 


and  (o  =  2;cf .  The  Polder  tensors  (31)  are  given  assuming  the  “engineer’s  time  convention”  of 

e+j«*  jhg  corresponding  tensors  for  the  “physicist’s”  time  convention  are  obtained  by  taking  the 
complex  conjugate  of  (31). 

In  (32)  and  (33),  several  different  important  frequencies  are  used.  The  operating  frequency  is  denoted 
by  f.  The  Lamor  or  precession  frequency  is  indicated  by  f0  while  the  forced  precession  frequency  is 
represented  by  fm .  The  precession  frequency  in  MHz  is  given  by 


f0MH2  =  2.8xH0  (34) 

4 

where  H0  is  the  applied  static  magnetic  bias  field  within  the  ferrite  sample  and  in  (34)  this  quantity  is 
assumed  to  be  given  in  terms  of  Oersteds  (Oe).  Note  that  H0  is  not  necessarily  equal  to  the  external 

applied  static  magnetic  field  (He)  since  for  practical  sample  geometries,  demagnitization  may  occur 
such  that  H0  <  He  depending  on  sample  geometry  and  applied  field  orientation.  The  interested 
reader  is  encourage  to  review  Pozar’s  text  book  chapter  concerning  ferrite  materials  for  further 
details1 .  The  forced  precession  frequency  in  MHz  is  given  by 

frz=2.8x(47tMs)  (35) 

where  the  saturation  magnetization  (4uMs)is  assumed  to  be  given  in  Gauss.  Msis  the  peak 
magnetization  achieved  within  a  ferrite  sample  as  an  increasingly  strong  static  bias  field  is  applied. 
Hence  it  represents  an  important  intrinsic  material  property  of  ferrites. 

The  co-pol  permeability  represented  by  (32)  and  (33)  includes  no  magnetic  loss.  To  include  such 
loss,  a  line  width  is  typically  used.  This  line  width  (AH)  is  illustrated  in  Figure  C-l 
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W,  N.  H2  tfo 

Figure  C-l .  Illustration  of  magnetic  line  width  (After  Pozar[l]). 
and  the  effect  is  to  modify  (34)  by  including  in  an  imaginary  component 

fJ“=Mx(H„+j^)  (36) 

Hence,  the  forced  precession  frequency  is  now  complex  rather  than  purely  real. 

An  important  parameter  for  total  electric  field  finite  element  analysis,  such  as  the  formulation 
presented  in  Section  B  is  the  inverse  permeability  tensor 

1  0  0 
0  A  'IVa  x-bias 

0  +iKA  %  . 

"  Vl  °  ~>A 

n;'  =  0  1  0  y-bias  (37) 

0  %  . 

Vl  -j k/a  o' 

+  j  /'a  %  0  z~bias 

0  0  1 

where  the  determinant  of  die  tensor  (3 1)  is  given  by  A  =  p.2  - k  2 . 

Several  interesting  phenomena  and/or  difficulties  arise  as  a  result  of  the  various  frequency  parameters 
associated  with  a  ferrite  material  (e.g.  f0 ,  fm ,  H0 ,  AH,  and  4tcM$ ).  For  example,  when  the  operating 
frequency  is  equal  to  the  precession  frequency  f  =  f0,  a  gyromagnetic  resonance  occurs  and  the 
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ferrite  permeability  is  infinite.  Since  the  skin  depth  is  inversely  proportional  to  the  permeability,  at 
the  gyromagnetic  resonance,  the  field  cannot  penetrate  into  the  ferrite  sample.  Alternatively,  the 
ferrite  sample  is  very  lossy.  This  condition  cannot  occur  when  a  magnetic  loss  is  included  since  the 
forced  precession  frequency  is  complex  while  the  operating  frequency  is  real.  Another  interesting 
phenomena  occurs  when  the  effective  permeability  (See  Section  E )  becomes  zero.  At  this  frequency, 
a  anti-ferromagnetic  resonance  (AFMR)  occurs  and  the  skin  depth  is  infinite.  Hence,  the  ferrite  is 
perfectly  transparent.  These  are  narrow  band  effects,  however,  there  is  some  hope  for  extending  these 
effects  over  a  wider  bandwidth  and  this  is  will  be  the  subject  of  future  research. 

A  difficulty  in  the  total  electric  field  finite  element  formulation  occurs  when  the  precession  frequency 
is  real.  At  some  frequency,  the  determinant  vanishes  (e.g.  when  |p|  =  |k|)  and  hence  the  inverse 

tensor  (37)  does  not  exist.  Once  again,  inclusion  of  a  small  magnetic  loss  (AH  *  0)  will  avoid  this 
problem.  In  addition,  as  UM  has  discovered,  the  material  model  presented  above  leads  to  solution 
difficulties  at  certain  frequencies  and  UM’s  findings  will  be  discussed  in  greater  detail  in  a  following 
section. 

Ferrite  materials,  due  to  their  Hermetian  (complex  conjugate,  transposed)  tensor  structure,  exhibit 
several  interesting  features.  The  polarization  properties  depend  on  the  orientation  of  the  applied  bias 
field  as  shown  in  (31)  and  (37).  For  example,  for  a  z-biased  condition,  ferrite  materials  “look” 
different  to  right-hand  circular  and  left-hand  circular  polarized  fields.  Hence,  the  propagation 
velocity,  loss,  and  other  physical  effects  are  different.  This  property  suggests  various  special  uses  for 
ferrites  in  antenna  design.  Alternatively,  depending  on  the  bias  state,  a  ferrite  material  can  look 
nearly  transparent  to  one  polarization  and  nearly  opaque  to  the  orthogonal  polarization. 

D.  Validation 

Both  MRC  and  UM  modified  their  respective  FE-BI  brick  computer  programs  to  accommodate 
anisotropic,  and  in  particular  ferrite,  materials.  Both  organizations  did  these  modifications  in  parallel 
to  assure  transfer  of  information  between  MRC  and  UM  and  to  provide  both  organizations  with  a 
first-hand  understanding  of  the  issues  involved  in  writing  FE-BI  computer  programs  for  ferrites.  This 
is  an  important  “learning  curve”  opportunity  for  both  MRC  and  UM  since  in  Phase  n,  both 
organizations  are  tasked  with  computer  program  development  that  will  eventually  be  merged  into  one 
common  program. 

The  specific  modifications  to  each  organization’s  software  were: 

1.  Implementing  the  anisotropic  finite  element  matrices  given  in  Section  B  (note,  the 
isotropic  material  matrices  are  obtained  by  setting  all  diagonal  terms  of  the  permittivity 
and  permeability  tensors  to  identical,  non-zero  numbers,  and  setting  all  other  tensor 
entries  to  zero). 

2.  Implementing  an  asymmetric  iterative  matrix  solver. 

3.  Automatic  calculation  of  the  permeability  tensor  for  ferrites  given  the  appropriate 
parameters. 


Mission  Research  Proprietary  Information 


Mission  Research  Proprietary  Information 


18 


MRC  implemented  a  biconjugate  gradient  (BiCG),  conjugate  gradient  squared  (CGS),  and 
biconjugate  gradient  stabilized  (BiCG-STAB)  iterative  matrix  solvers.  The  BiCG  method  was  found 
to  be  the  preferred  algorithm  in  terms  of  convergence  rate  and  ability  to  reach  a  converged  condition. 
The  solution  generated  using  the  BiCG  algorithm  and  a  LU  decomposition  solver  (UNPACK)  is 
shown  in  Figure  D-l . 


Figure  D-l .  Agreement  between  BiCG  and  LU  decomposition  solvers. 

Although  this  particular  geometry  exhibited  solution  difficulties  near  the  forced  precession  frequency 
(2.4  GHz  for  this  example),  the  agreement  between  the  BiCG  and  LU  decomposition  solvers 
demonstrates  the  validity  of  the  BiCG  solver  as  implemented  by  MRC  for  asymmetric  matrices. 

The  next  validation  task  was  to  compare  MRC’s  and  UM’s  computer  programs  with  published 
reference  data.  This  was  a  difficult  task  for  two  reasons: 

1.  The  published  data  rarely  corresponded  to  a  cavity-backed  geometry  (which  is  what  the 
FE-BI  computer  program  assumed)  and  hence  any  comparisons  involved  potentially 
important  geometrical  approximations. 

2.  Published  solutions  often  were  not  consistent  with  one  another,  likely  due  to  the  special 
difficulties  in  modeling  ferrite  material. 

Various  sources  were  used  and  the  comparisons  along  with  a  brief  discussion  of  each  will  be 
presented  in  turn. 
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D.l  Cavity-backed  Apertures 

One  source  of  reference  data  that  was  particularly  useful  was  found  in  Dr.  David  KokotofFs  Ph.D. 
dissertation2.  In  this  work.  Dr.  Kokotoff  developed  a  method  of  moments  (MoM)  computer  program 
for  cavity-backed,  ferrite  loaded  apertures.  This  is  essentially  the  MoM  equivalent  of  MRC’s  FE-BI 
computer  program,  except  the  MoM  program  is  only  capable  of  modeling  layered  cavity  fills  as 
opposed  to  fully  inhomogeneous  fills.  Nevertheless,  the  data  published  in  [2]  was  useful  for 
comparison.  Figure  D-2  illustrates  a  geometry  used  by  both  Kokotoff  and  MRC. 


Figure  D-2.  Layered  cavity-backed  aperture  used  for  comparison  between  MoM  and  FE-BI. 

Kokotoff  ran  his  MoM  computer  program  for  various  bias  states  (400  Oe  <  H0  <  700  Oe)  and 
recorded  the  radar  cross  section  (RCS)  within  a  300  MHz  bandwidth.  IBs  results  are  reproduced  in 
Figure  D-3. 


Figure  D-3.  Reproduction  of  Figure  6. 1  in  [2]. 

MRC’s  results  for  the  same  geometry  (see  Figure  D-2)  are  presented  in  Figure  D-4. 
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Frequency  [G Hz] 


Figure  D-4.  MRC's  FE-BI  comparison  to  Figure  D-3. 

The  agreement  is  quite  good.  This  result  is  significant  in  that  it  demonstrates  the  validity  of  the  FE- 
BI  implementation  by  comparison  with  a  totally  different  formulation. 

Not  all  comparisons  between  KokotofFs  results  and  either  MRC’s  or  UM’s  results  were  favorable. 
Tremendous  difficulty  was  encountered  for  certain  frequency  ranges  and  both  MRC  and  UM 
conclude  that  in  these  bands,  the  ferrite  material  exhibits  behavior  that  is  ill-favorable  for  electric  field 
FE-BI  analysis.  For  example,  in  the  vicinity  of  the  forced  precession  frequency  (and  an  unbiased 
ferrite  sample)  as  mentioned  previously,  the  determinant  of  the  permeability  tensor  goes  to  zero  and 
hence  the  inverse  permeability  tensor  is  not  defined.  This  results  in  erratic  behavior  and  ill- 
conditioned  matrices.  Figure  D-5  illustrates  the  condition  number  for  a  dielectric-  and  ferrite-filled 
cavity  problem.  In  this,  the  dielectric  (isotropic)  condition  number  is  generally  good  and  even  its 
peak  near  2.2  GHz  is  not  too  bad.  However,  clearly  the  ferrite-filled  case  exhibits  poor  conditioning 
just  below  the  aperture  resonance  (~2.2  GHz)  and  near  the  forced  precession  frequency  (~2.4  GHz). 
For  this  case,  the  bias  field  was  shut  off  and  no  magnetic  loss  was  assumed,  so  the  determinant  of  the 
permeability  tensor  vanished  at  2.4  GHz. 
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Frequency  [GHz] 

Figure  D-5.  Condition  number  verses  frequency  (large  condition  numbers  indicate  a  poorly 
conditioned  matrix). 

D.2  Patch  Antennas 

Patch  antennas  are  another  class  of  antennas  investigated  under  this  effort.  Unfortunately,  the 
available  reference  data  for  ferrite-backed  patch  antennas  assume  an  infinite  ferrite  substrate  whereas 
the  FE-BI  computer  program  used  by  MRC  and  UM  assume  a  finite  cavity  below  the  patch.  Figure 
D-6  reproduces  a  RCS  calculation  for  a  ferrite-backed  patch  antenna  presented  by  Pozar3. 


2.0  2.3  3.0  S.S  4.0  4.9  S.O  S.S  3.0  3.5  7.0  7.3  3.0 

Fraquancy  (GHz) 

Figure  D-6.  Reproduction  of  Pozar’s  [3]  results. 
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Figure  D-7  illustrates  MRC’s  results  for  the  same  patch  antenna  (except  the  substrate  is  recessed  in  a 
cavity  rather  than  being  an  infinite  grounded  ferrite  slab). 


Figure  D-7.  Scattering  by  a  1.3xl.3  cm  patch  antenna  recessed  in  a  2.6x2.6x0.13  cm  ferrite  filled 
cavity. 

The  agreement  is  good  except  for  low  frequencies  where  the  cavity-backed  structure  is  not  a  good 
approximation  to  the  original  grounded  ferrite  slab.  This  is  due  to  the  fact  that  below  resonance, 
aperture  edges  are  not  electrically  far  from  the  patch  and  hence  their  presence  has  a  greater  effect  on 
the  scattered  fields  attributed  to  the  patch  antenna. 

E.  Convergence  Issues 

Several  factors  affect  the  convergence  properties  of  iterative  matrix  solvers.  These  include: 

1.  Aspect  ratio  of  elements  (e.g.  large  differences  in  element  shapes  throughout  the  mesh  is 
bad). 

2.  Excitation  (typically  radiation  analysis  requires  more  iterations  than  RCS  analysis). 

3 .  Formulation  (FE-BI  is  generally  better  conditioned  than  FE-ABC,  FE-AA,  or  FE-PML). 

4.  Material  properties  (anisotropic  is  often  more  difficult  to  solve  than  isotropic  materials). 

This  section  concentrates  on  the  last  item  since  it  is  particular  to  ferrite  materials*. 

E.1  Propagation  Constants 

Although  convergence  was  achieved  for  the  previous  examples,  convergence  difficulties  were 
encountered  while  studying  the  ferrite  materials.  To  examine  the  reason  behind  these  convergence 
difficulties,  it  is  necessary  to  look  at  the  wave  modes  which  are  excited  within  the  ferrite  region. 
These  modes  are  linked  to  the  propagation  constants  within  the  ferrite  material.  A  correlation  can  be 
made  between  the  propagation  constants  and  the  regions  of  non-convergence. 


1  This  section  is  taken  verbatim  from  a  report  submitted  by  UM  to  MRC  in  April  1997. 
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When  dealing  with  ferrite  materials,  the  field  behavior  is  determined  by  the  propagation  direction  and 
its  orientation  with  the  applied  magnetic  bias  field  direction.  There  are  two  separate  cases  which 

determine  the  effective  permeability  (p^)  within  the  ferrite  -  the  longitudinal  case  where  the 
propagation  is  parallel  to  the  applied  bias  field  and  the  transverse  case  where  propagation  is 
perpendicular  to  the  applied  bias  field.  In  the  longitudinal  case 


peff  =  p±K  (38) 


whereas  in  the  transverse  case 


Heff  = 


2  2 

p  -k 


ft 


(39) 


For  both  instances,  the  propagation  constant  within  the  ferrite  is  calculated  as 


Y  =  j©7eoPoerH«fr  =jk0  Ap^  =a  +  jP  (40) 

t 

A  careful  study  of  die  solver  convergence  rate  for  various  propagation  constants  was  done.  It 
provided  the  conclusion  that  convergence  behavior  can  be  predicted,  by  looking  at  each  of  the  ferrite 
modes  and  the  associated  p  and  k  parameters. 


E.1.1  Ferrite  Cavity  Example 

To  demonstrate  die  correlation  between  y  and  convergence,  the  geometry  shown  in  Figure  E-l  was 
considered. 


X 


y 


Layer!  d  =  1.905  cm,  £r  =  1.0,  p*  -  1.0 

Layer?  d  =  0.635  cm.  &  =  13.9.  4kM*=BOOG.  HosIOOOo,  AH=1  OO© 

Layer3  d  =  1 .905  cm,  &  =  1 .0,  jir  =  1 J) 

Figure  E-l.  Ferrite  cavity  geometry. 

In  this  example,  H0  =  100Oeand  Figure  E-2  and  Figure  E-3  show  the  behavior  of  pand  k, 
respectively,  verses  frequency. 
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FigureE-2.  pvs.  f,  H0  =1000e,  4teMs  =800G,  AH  =  1 0 Oe . 


We  observe  that  p  is  less  than  unity  for  all  frequencies  in  the  desired  range,  and  this  situation  is 
similar  to  the  dual  case  of  sr  less  than  one.  The  latter  is  associated  with  plasma  behavior  of  the 
ionosphere.  At  those  frequencies  where  er  is  less  than  one,  the  ionosphere  becomes  completely 
reflective  and  gives  rise  to  the  phenomenon  of  short  wave  propagation  around  the  earth. 

Figure  E-4  and  Figure  E-5  show  a  and  P  for  three  different  values  of  p^ .  Excluding  the  mode 
with  p^  =  p  —  k  ,  the  other  two  modes  have  P  nearly  zero  for  all  frequencies.  More  specifically  P 
exhibits  a  filter-like  behavior  from  1.3  -  2.5  GHz  as  demonstrated  in  Figure  E-5  for  the  extraordinary 
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(transverse)  mode.  Interestingly,  the  convergence  curve  in  Figure  E-6  shows  that  the  convergence  is 
very  slow  or  not  achieved  in  the  band  where  p  ~  0. 


|i  -K 

Figure  E-4.  a  normalized  to  k0,  -  |aeff  = - =  |^+K,--peff  =  |i-k 


\i  -K 

Figure  E-5.  p  normalized  to  k0,  -|reff  = - , —  |ieff  =^+K,--peff  =  (J.-K 
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Figure  E-6.  Convergence  performance,  H0  =  100  Oe,  #  of  unknowns  =  3915. 

E.1.2  Ferrite  Cavity  with  Patch  Antenna  Example 

Figure  E-7  and  Figure  E-8  show  the  frequency  behavior  of  f x  and  k  .  From  the  figures,  both  p  and 
k  can  have  negative  values  for  the  frequency  range  considered.  Also  Figure  E-9  and  Figure  E-10 
show  that  the  attenuation  and  phase  constants  exhibit  a  behavior  similar  to  the  ferrite  cavity.  Once 
again,  the  convergence  performance  (See  Figure  E-l  1)  is  governed  by  the  behavior  of  the  transverse 
P .  When  P  is  nearly  zero,  we  again  see  that  no  convergence  is  achieved.  More  specifically,  the 
frequencies  at  which  the  transverse  P  is  zero  corresponds  to  the  frequency  band  of  non-convergence. 


Figure  E-7.  pvs.  f,  H0  =6000e,  47tMs  =650G,  AH=30Oe. 
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F.  Phase  I  Trade  Studies 

MRC  and  UM  conducted  a  number  of  trade  studies  during  Phase  I  to  establish  the  feasibility  of 
applying  finite  element  analysis  to  ferrite  materials  and  to  aid  in  the  definition  of  a  Phase  II 
development  plan.  This  section  presents  the  results  of  these  trade  studies. 

F.l  Perfectly  Match  Layer  Mesh  Truncation 

The  University  of  Michigan  has  been  studying  the  use  of  perfectly  matched  layer  (PML)  mesh 
truncation  condition.  Their  work  has  involved  the  anisotropic  PML  formulation  first  presented  by 
Sacks,  et  al.4.  Figure  F-l  illustrates  this  anisotropic  PML.  Use  of  this  PML  on  curved  boundaries  has 
not  been  established  and  therefore  was  investigated  by  UM.  In  addition,  since  this  PML  has  various 
parameters,  a  trade  study  to  determine  the  optimal  parameter  mix  was  undertaken.  Figure  F-2 
illustrates  the  results  of  this  trade  study.  From  this  chart,  an  absorber  designer  can  determine  the 
optimal  parameters  which  result  in  a  minimally  thick  absorber  with  maximum  absorption  for  a  given 
number  of  absorber  layers. 


Figure  F-l .  Comparison  of  new  anisotropic  PML  and  traditional  PML. 


Mission  Research  Proprietary  Information 


Mission  Research  Proprietary  Information 


30 


Figure  F-2.  PML  absorber  design  curves  where  0t/>,o  is  the  thickness  of  the  PML  and  N  is  a  PML 
parameter  describing  the  material  properties.  The  straight  lines  give  |R|  in  dB  and  the  curved  lines 

give  N:  ( — )  exact  |R|,  (•  •  •)  homogeneous  case  and  ( - )  inhomogeneous  case. 


F.1.1  Curved  PML  Trade  Study 

UM  undertook  a  trade  study  to  determine  the  performance  of  the  PML  for  curved  surfaces.  Their 
overall  conclusion  is  that  this  PML  will  work  well  for  a  curved  surface;  however,  not  as  well  as  in  the 
case  of  a  flat  surface.  This  PML  is  superior  to  more  traditional  second  order  absorbing  boundary 
conditions  (ABCs).  Figure  F-3  illustrates  UM’s  results. 
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Figure  6:  S  for  ABCs  and  material  absorbers  with  xo  =  10,  d  =*  Ao/2  and  it  =*  1  —  j'kj: 

( - )  ABCs,  ( . )  homogeneous  isotropic  material,  ( - •)  homogeneous  anisotropic 

and  (—  — )  inhomogeneous  anisotropic  material.  For  the  materials,  the  top  curve 
is  for  a  layer  with  r  —  0.18Aq  and  the  bottom  one  is  for  an  infinite  medium. 


Figure  7:  6  for  analytical  and  FEM  results  with  xo  =  10,  d  —  Ao/2,  r  —  O.ISAq  and 
«  » 1  -  jiq:  ( - )  analytical  and  (♦  ♦  ♦)  FEM. 

Figure  F-3.  Results  of  the  curved  PML  trade  study. 

F.2  Matrix  Iterative  Solvers  and  Conditioning 

UM  undertook  a  study  of  various  iterative  matrix  solvers  due  to  MRC’s  preliminary  discovery  that 
ferrite  materials  have  a  rather  detrimental  effect  on  the  matrix  condition  number.  UM  found  that  the 
GMRES  solver  required  fewer  iteration  than  either  the  BiCG  or  QMR  algorithms  (this  is  to  be 
expected  based  upon  published  results  in  the  literature).  Figure  F-4  illustrates  the  convergence 
behavior  of  these  three  iterative  solvers  (QMR,  BiCG,  and  GMRES)  for  a  3000  unknown  problem. 
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atptm«4>eta,  Dlag  Precond. 


Figure  F-4.  Convergence  behavior  for  three  iterative  solvers  (3000  unknowns). 

A  study  of  the  benefits  of  using  matrix  preconditioning  was  also  undertaken.  Figure  F-5  illustrates 
the  benefits  of  such  preconditioning.  Preliminary  indication  from  both  MRC  and  UM  research  staff 
indicates  that  preconditioning  is  essential  to  achieving  reliable  and  accurate  result  from  a  finite 
element  analysis  when  ferrite  materials  are  present. 


BCG  solver,  alpha=beta 


50  1  00  1  50  2  0  0  2  5  0  30  0  3  5  0  400  4  5  0  500 

Iterations 


Figure  F-5.  Convergence  behavior  of  an  iterative  solution  with  and  without  preconditioning. 
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F.3  Element  Distortion 

UM  also  investigated  the  effect  of  element  distortion  on  condition  number.  It  is  a  well  known  fact 
that  the  best  conditioned  matrices  result  from  elements  that  have  roughly  a  1 : 1 : 1  aspect  ratio.  That  is, 
all  edges  of  the  element  are  approximately  equal.  Often  it  is  not  possible  to  obtain  such  elements 
without  a  high  computational  cost  in  terms  of  number  of  unknowns  (e.g.  lots  of  small  elements). 
Distortion  of  the  elements  is  expected  to  be  especially  important  for  ferrite  materials  due  to  their 
anisotropic  permeability.  Figure  F-6  illustrates  the  cases  investigated  by  UM  and  Figure  F-7  presents 
UM’s  results. 


Figure  F-6.  Element  distortion  cases  examined  (data  shown  in  Figure  F-7). 
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oondition  number  in  power  of  system  size 


condition  number  normalized  to  isotropy/ no rvd istortion 


Figure  F-7.  Effect  of  element  distortion  on  matrix  condition  number.  The  cases  correspond  to  the 
four  brick  images  shown  in  Figure  F-6. 

G.  Ferrite-backed  Patch  Antenna  Examples 

MRC  and  UM  exercised  the  FE-BI  computer  programs  developed  during  Phase  I  to  investigate  the 
radiation  properties  of  ferrite-backed  antennas.  This  section  presents  a  sample  of  our  results. 
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G.l  Radiation  Properties  of  a  Single  Ferrite-backed  Patch  Antenna 

As  a  partial  attempt  to  answer  the  question,  “What  good  are  these  antennas?”,  MRC  has  investigated 
the  radiation  properties  of  a  particular  patch  antenna  using  our  FE-BI  computer  program, 
LM  BRICK.  The  antenna  chosen  for  initial  study  was  the  one  used  by  Prof.  Pozar  to  demonstrate 
polarization  diversity  [3],  This  antenna  was  0.61  x  0.61  cm  printed  on  a  0.127  cm  thick  ferrite 
substrate.  The  substrate  had  a  dielectric  constant  of  15  and  a  magnetization  saturation  of  650  Gauss. 
The  substrate  was  normally  (z-axis)  biased. 

For  our  study,  we  placed  the  antenna  cited  above  in  a  1.22  x  1.22  cm  metallic  cavity  and  considered 
the  transverse  (x-  and  y-axis)  as  well  as  normally  biased  (z-axis)  cases.  Pozar’ s  integral  equation 
formulation  allowed  investigation  of  only  the  normally  biased  case.  We  found  that  all  three  bias 
orientation  have  useful  properties.  We  assumed  a  feed  point  at  (0,-0.3)  cm  relative  to  the  center  of  the 
patch.  This  was  close  to  an  edge  and  thereby  resulted  in  a  rather  high  resonant  resistance.  This  feed 
point  can  be  moved  closer  to  the  center  of  die  patch  to  reduce  the  resonant  resistance  for  matching 
purposes. 

We  varied  the  bias  field  strength  to  verify  the  resonant  frequency  tuning  reported  in  Pozar’ s  paper. 
Our  results  agreed  with  the  reported  data  (specifically.  Figure  2  of  Pozar’ s  paper  which  is  reproduced 
below  as  Figure  G-l)  within  a  small  frequency  shift.  Typically  such  shifts  can  be  reduced  by 
increased  sampling  of  the  geometry  volume;  however,  for  our  purposes,  our  results  agreed  with 
Pozar’s  to  a  sufficient  degree. 

For  the  remainder  of  this  section,  we  assume  a  bias  field  strength  of 400  Oe. 


Figure  G-l.  Figure  2  from  Pozar's  paper  [3],  This  illustrates  the  variation  of  bias  frequency  as  the 
bias  field  strength  is  changed  for  a  normally  biased  ferrite-backed  patch  antenna. 
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G.1.1  Dielectric-backed  Patch  Antenna 

MRC  initially  investigated  the  radiation  properties  of  a  dielectric-backed  patch  antenna  with  the  same 
size  as  Pozar’s  antenna  cited  above  and  with  a  dielectric  constant  of  15.  Figure  G-2  illustrates  the 
impedance  spectrum  of  such  and  antenna. 


0.61x0.61  Patch  Antenna  [Dielectric,  er=15] 


Frequency  [GHz] 


Figure  G-2.  Impedance  spectrum  of  a  0.61x0.61  cm  dielectric-backed  patch  antenna. 


In  Figure  G-2,  the  dielectric-backed  antenna  is  linearly  polarized  in  the  y-direction.  Had  we  fed  this 
antenna  at  (-0.3, 0.0)  cm  rather  than  (0.0, -0.3),  we  would  expect  similar  results;  however,  the  radiation 
would  now  be  polarized  in  the  x-direction. 

G.1.2  X-axis  Ferrite  Biasing 

MRC  modeled  a  ferrite-backed  antenna  where  the  bias  orientation  was  along  the  x-axis  (the  bias  field 
was  400  Oe).  Figure  G-3  illustrates  the  impedance  spectrum  for  this  antenna. 
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0.61x0.61  Patch  Antenna  [H0=4OO,4nMs=65O,x-axis,Er=15] 


Figure  G-3.  Impedance  spectrum  for  a  0.61x0.61  cm  patch  printed  on  a  x-axis  biased  ferrite. 

In  this  case,  the  resonant  frequency  was  5.65  GHz.  Note  that  it  resonates  at  a  slightly  higher 
frequency  due  to  the  application  of  a  bias  field.  For  a  transversely  biased  ferrite,  the  polarization  is 
linear  and  in  this  case  the  antenna  is  y-polarized. 

G.1.3  Y-axis  Ferrite  Biasing 

MRC  modeled  a  ferrite-backed  antenna  where  the  bias  orientation  was  along  the  y-axis  (the  bias  field 
was  400  Oe).  Figure  G-4  illustrates  the  impedance  spectrum  for  this  antenna.  In  this  case,  the 
resonant  frequency  is  6.05  GHz.  Notice  the  rather  large  shift  in  resonant  frequency  for  this  transverse 
bias  orientation  as  compared  to  the  dielectric  (5.61  GHz)  and  the  x-axis  bias  case  (5.65  GHz).  The 
radiation  pattern  is  again  linear  and  the  radiated  field  is  polarized  in  the  y-direction.  The  cross- 
polarization  term  is  asymmetric  due  to  the  position  of  the  feed. 

Note  that  for  transversely  biased  antennas  (e.g.  x-  or  y-axis),  the  polarization  is  linear,  however,  the 
resonant  frequency  can  be  tuned  over  two  ranges  by  varying  the  bias  field.  The  two  ranges  are 
selected  by  the  bias  field  orientation,  hi  the  cases  presented  herein,  y-directed  polarization  was 
obtained  by  feeding  the  antenna  along  the  y-axis.  If  we  desire  x-directed  polarization,  placing  the 
feed  on  the  x-axis  would  be  appropriate. 
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0.61x0.61  Patch  Antenna  [H0=4OO,4nMs=65O,y-axis,Er=15] 


Frequency  [GHz] 


Figure  G-4.  Impedance  spectrum  of  a  0.61x0.61  cm  patch  antenna  on  a  y-biased  ferrite. 

G.1.4  Z-axis  Ferrite  Biasing 

MRC  modeled  a  ferrite-backed  antenna  where  the  bias  orientation  was  along  the  z-axis  (the  bias  field 
was  400  Oe).  Figure  G-5  illustrates  the  impedance  spectrum  for  this  antenna.  In  this  case,  two 
resonant  frequencies  are  observed  (as  reported  by  Pozar  and  shown  in  Figure  1).  One  resonant 
frequency  is  at  5.1 15  GHz  while  the  other  is  at  6.445  GHz. 


0.61x0.61  cm  Patch  Antenna  [H0=4OO,4icMs=65O,z-axis,£r=15] 


Figure  G-5.  Impedance  spectrum  for  a  0.61x0.61  cm  antenna  on  a  normally  biased  ferrite. 
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The  radiation  pattern  at  5.115  GHz  is  left-hand  circular  polarized  (LHCP)  and  it  is  shown  for  the  x- 
and  y-axis  cuts  in  Figure  G-6.  The  corresponding  pattern  at  6.445  GHz,  which  is  right-hand  circular 
polarized  (RHCP),  is  shown  in  Figure  G-7. 

0.61x0.61  cm  Patch  Antenna  [H0=4OO,4JtMs=65O,z-axis,ef=15] 


Angle  [©.deg] 


Figure  G-6.  Gain  pattern  for  a  normally  biased  patch  antenna  at  5. 1 1 5  GHz. 
0.61x0.61  cm  Patch  Antenna  [H0=4OO,4JlMs=65O,z-axis,Er=15] 


Angle  [0,deg] 

Figure  G-7.  Gain  pattern  for  a  normally  biased  patch  antenna  at  6.445  GHz. 
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G.1.5  Summary  of  Radiation  Simulations 

Clearly,  a  normally  biased  antenna  can  be  operated  in  one  of  two  circular  polarization  modes  at  two 
distinct  frequencies.  Also,  linear  polarization  is  obtained  when  the  bias  field  is  transverse  to  the 
ferrite  (again  at  two  distinct  frequencies).  Fine  frequency  tuning  can  be  accomplished  by  varying  the 
bias  field  strength. 

These  results  suggest  that  a  polarization  and  frequency  diverse  square  antenna  element  can  be 
obtained  by  providing  two  feed  points  (one  along  the  x-axis,  one  along  the  y-axis)  only  one  of  which 
is  excited  for  a  desired  polarization.  A  means  of  the  bias  field  strength  and  orientation  must  be 
provided  perhaps  through  the  use  of  small  electromagnets.  This  antenna,  with  no  moving  parts, 
would  be  capable  of  various  linear  polarization  states  as  well  as  a  circular  polarization  mode.  It 
would  also  be  capable  of  operating  over  a  large  bandwidth  by  altering  the  bias  field  direction  and 
strength  even  though  the  instantaneous  bandwidth  is  quite  small.  MRC  is  currently  investigating  the 
use  of  two  feeds  at  the  same  time;  however,  the  feeds  are  fed  90  degrees  out-of-phase.  We  anticipate 
circular  polarization  for  the  dielectric  and  z-axis  biased  ferrite  while  the  x-  and  y-axis  biased  ferrite 
configurations  will  remain  linearly  polarized.  This  is  due  to  the  fact  that  the  two  feed  points  in  the 
case  of  x-  or  y-biased  ferrite  patch  antennas  excite  different  antenna  resonances  and  hence  at  most 
only  one  feed  or  the  other  will  be  “active”  (e.g.  matched)  and  any  given  frequency. 

G.2  Large  Ferrite-backed  Antenna  Array 

MRC  modeled  a  1 1  x  1 1  patch  antenna  array  on  a  Silicon  Graphics  workstation  using  only  60  MB  of 
RAM  and  double  precision  calculations.  Each  antenna  element  was  4  x  3  cm  and  placed  on  top  of  a  8 
x  6  x  0. 15  cm  ferrite  slab  (within  a  metal  cavity).  The  center-to-center  element  spacing  was  10  cm  in 
both  directions  yielding  a  total  array  size  of  approximately  lxl  meter  and  resulting  number  of 
unknowns  was  approximately  50,000.  The  gain  pattern  for  this  antenna  is  shown  in  Figure  G-8.  In 
this,  the  main  lobe  was  steered  to  30  degrees  off  broadside  using  standard  beam  steering  techniques. 
Computation  of  this  pattern  required  approximately  8.5  hours  on  a  Silicon  Graphics  workstation. 

To  model  a  similar  antenna  using  a  second-order  ABC  and  brick  elements,  MRC  estimates 
approximately  775,000  unknowns  (ABC  is  five  layers  from  antenna  +  50,000  unknowns  for  the 
antenna  itself).  Clearly,  such  an  explosive  growth  in  the  number  of  unknowns  is  not  conducive  to 
accurate  and  efficient  solution.  Assuming  an  average  bandwidth  of  60,  a  direct  solver  used  in  a  FE- 
ABC  computer  program  would  require  approximately  744  MB  of  RAM  before  fill-in!  We  note  that 
for  similar  zoning,  the  number  of  unknowns  in  a  tetrahedral  mesh  will  be  higher  than  the  bride  mesh 
cited  above. 
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Angle  [deg] 

Figure  G-8.  Radiation  pattern  from  an  11  x  11  array  of  patch  antennas  printed  on  a  ferrite  substrate. 
This  example  illustrates  MRC's  ability  to  model  large  structures  on  ferrite  materials  including 
antenna  elements  and phase  shifting  devices. 

H.  Non-Uniform  Static  Bias  Field  Effects 

Practical  ferrite  antenna  construction  often  involves  a  permanent  magnet  to  provide  die  required  static 
magnetic  bias  field.  Most  CEM  analysis  techniques  for  ferrite  material  (in  particular  the  method  of 
moments)  assume  implicitly  a  uniform  bias  field  distribution.  However,  since  permanent  magnets  are 
finite,  the  static  field  attributed  to  them  is  non-uniform.  Figure  H-l  illustrates  the  static  magnetic 
field  within  a  metallic  cavity  as  measured  by  Dr.  Kokotoff  [2],  Note  that  the  field  has  roughly  a 
parabolic  shape  and  that  introduction  of  steel  bars  between  the  magnet  and  the  ferrite-filled  cavity 
reduces  the  distortion  in  the  static  field.  Also,  as  the  static  field  increases  (here  by  use  of  more 
permanent  magnets),  the  “bowing”  in  the  static  field  becomes  more  distinct. 

The  non-uniform  field  encountered  in  practice  encourages  the  use  of  the  finite  element  method  since 
the  non-uniform  bias  field  results  in  an  magnetically  inhomogeneous  ferrite  sample. 
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Figure  H-l .  Measured  magnetic  field  distribution  within  cavity  from  [2], 

As  an  example  of  the  effect  of  a  non-uniform  bias  field,  consider  the  scattering  bya6cmx6cmxl 
cm  ferrite  filled  cavity.  Figure  H-2  illustrates  the  variation  in  RCS  as  a  function  of  frequency  and 
non-uniform  bias  state.  Note  that  the  RCS  near  1.5  GHz  is  prominently  affected  by  the  bias  field 
distribution  throughout  the  cavity. 

RCS  of  a  6x6x1  cm  Cavity  w/ Various  Biasing  Schemes 


Frequency  [G Hz] 

Figure  H-2.  Illustration  of  non-uniform  bias  field  effect  on  the  RCS  of  a  ferrite-filled  cavity. 
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This  example  demonstrates  the  requirement  for  a  finite  element  computer  program  since  MoM  cannot 
model  non-uniform  bias  fields.  Further  examples  of  the  importance  of  including  non-uniform  bias 
fields  will  be  presented  in  Section  I. 

I.  Ferrite  Radomes:  An  Electronic  Shutter  Concept 

An  important  proposed  application  of  ferrite  materials  is  as  a  radome.  This  is  due  to  the  fact  that  the 
magnetically  tunable  properties  (loss  and  effective  permeability)  are  potentially  useful  as  an  electronic 
shutter.  Typically,  radomes  are  designed  to  permit  radiation  from  the  antenna  to  the  outside  world 
either  at  all  frequencies  (a  dielectric  radome)  or  in  a  limited  band  of  frequencies  using  a  frequency 
selective  surface  (FSS).  Unfortunately,  the  reality  of  physics  is  that  a  good  antenna  is  a  good 
scattering  body  and  hence  if  a  radome  is  transparent,  the  RCS  of  the  antenna  can  be  a  significant 
factor  in  the  over-all  platform  RCS.  Hence  a  need  to  shutter  the  radome  and  make  it  transparent  only 
during  active  radiation. 

MRC  undertook  at  the  end  of  Phase  I  a  brief  feasibility  study  concerning  the  use  of  ferrites  as  an 
electronic  shutter.  The  results  are  somewhat  promising.  The  ferrite  layer  does  indeed  act  like  a 
shutter  radome,  though  it  is  very  narrow  banded.  A  solution  for  this  narrow  bandwidth  of  operation 
should  become  an  active  research  topic  in  order  to  utilize  ferrites  as  shutter  radomes.  This  section 
presents  MRC’s  investigation  of  ferrite  radomes  beginning  with  the  modified  FE-BI  computer 
program  necessary  for  determining  the  transmission  coefficient  of  a  ferrite  sample. 

LI  Two  Aperture  FE-BI  Computer  Program  for  Transmission  Calculations 

MRC  modified  the  FE-BI  computer  program  discussed  in  Section  B.  This  modification  consisted  of 
opening  an  aperture  at  the  base  of  the  metallic  cavity  and  hence  providing  the  capability  of  modeling 
thick  metallic  planes.  The  boundary  integral  submatrix  could  be  re-used  since  the  discretization  is 
identical  for  the  upper  and  lower  apertures  when  brick  elements  are  used  in  the  volume.  This 
computer  program  requires  a  few  more  unknowns  but  no  new  calculation  or  storage  of  the  matrix 
entries  hence  its  efficiency  is  comparable  to  the  single  aperture  brick  FE-BI  computer  program. 

L2  Transmission  Coefficient  Calculations 

The  transmission  coefficient  for  a  ferrite  sample  was  calculated  by  first  computing  the  far-zone 
forward  scattered  fields  when  the  volume  between  the  apertures  is  filled  with  air.  Then  the  forward 
scattered  fields  are  computed  with  the  ferrite  material  in  the  cavity.  The  transmission  coefficient  was 
then  calculated  by  dividing  the  latter  by  the  former 


where  the  fields  are  complex  and  the  division  is  done  on  a  component-by-component  basis.  This 
scheme  removes  the  fixture  scattering  (e.g.  thick  metal  plane  with  an  aperture)  from  the  calculation  to 
isolate  the  ferrite  effects.  Figure  1-1  illustrates  the  transmission  coefficient  for  a  ferrite  sample  with 
various  uniform  and  non-uniform  bias  field  distributions. 
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Non-uniform  Bias  Field  Effect  [4nMs=1780  G] 


Frequency  [GHz] 

Figure  1-1 .  Transmission  coefficient  for  a  ferrite  substrate  under  various  bias  conditions. 

Clearly,  the  transmission  coefficient  is  determined  by  the  bias  field  and  a  non-uniform  bias  field  has 
an  effect.  The  transmission  null  (e.g.  when  the  radome  is  opaque)  is  a  characteristic  of  ferrites  and  its 
location  is  determined  by  the  ferrite  material  properties  and  bias  field  distribution. 

J.  Concluding  Remarks  and  Suggestions  for  Future  Work 

jfTxning  Phase  I  of  this  SBIR  effort,  MRC  and  UM  demonstrated  the  feasibility  of  using  modem  finite 
element  analysis  to  characterize  the  behavior  of  ferrite  materials.  Specifically,  MRC  and  UM 
investigated  the  following: 

1 .  Finite  element  formulations  for  ferrite  materials. 

2.  Implementation  of  those  formulas  in  a  computer  program. 

3 .  Convergence  issues  when  ferrite  materials  are  being  modeled. 

4.  Mesh  truncation  schemes. 

5.  Element  distortion  effects. 

6.  and  uses  of  ferrite  antennas  as  frequency-agile  antennas,  polarization-agile  antennas  and 
electronic  shutter  radomes.'| 

The  importance  of  these  studies  will  be  realized  in  Phase  II  as  the  MRC  Team  develops  a  state-of-the- 
art  conformal  ferrite  antenna  modeling  package. 

t 

J.l  Suggestions  for  Future  Work 

The  MRC  Team  suggests  the  following  engineering  services  relating  to  finite  element  analysis 
research  and  development  be  performed: 

1.  Modify  existing  Phase  I  brick  finite  element  program  to  incorporated  advances  made 
during  Phase  II  in  terms  of  iterative  equation  solvers  or  optimization  techniques. 
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2.  Develop  a  tetrahedral  finite  element  program  capable  of  modeling  ferrites,  dielectrics,  and 
perfect  electric  conductors. 

3.  Develop  models  of  the  static  magnetic  field  within  the  antenna  cavity  to  include  non- 
uniform  bias  field  effects. 

4.  Develop  user-oriented  antenna  design  tools  that  incorporate  advanced  frequency 
extrapolation  and  antenna  optimization  features. 

5.  Investigate  improved  matrix  solution  algorithms. 

6.  Investigate  improved  (more  flexible)  mesh  truncation  schemes. 

7.  Investigate  ferrite  antenna  uses  through  a  combination  of  modeling,  prototype  fabrication, 
and  measurement. 

8.  Investigate  novel  ferrite-based  feed  networks  that  will  offer  advantages  in  toms  of  size, 
weight,  and  reduced  complexity  as  compared  to  current  phased  array  feed  networks. 
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